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^ I The response to the field sequence of nonresonant hole burning, a pump-wait-probe experi- 

^ I ment originally designed to investigate slow relaxation in complex systems, is calculated for 

CN I a model of Brownian oscillators, thus including inertial effects. In the overdamped regime 

bJOl the model predictions are very similar to those of the purely dissipative stochastic models 

investigated earlier, including the possibility to discriminate between dynamic homogeneous 
and heterogeneous relaxation. The case of underdamped oscillations is of particular interest 
when low-frequency excitations in glassy systems are considered. We show that also in this 
^ situation a frequency selective modification of the response should be feasable. This means 

\^ • that it is possible to specifically address various parts of the spectrum. An experimental 

in 
o 

00 ■ further light on the nature of the vibrations around the so-called boson peak. 
O 

cn , 

O . PACS Numbers: 64.70 Pf,05.40.+j,61.20.Lc 

a 

^ . 

§ ■ I. Introduction 
o 

The relaxation functions observed in disordered materials such as glasses, spin-glasses, 
^ i disordered crystals or proteins are usually found to decay non-exponentially on macroscopic 
^ I time scalesPP. Several experimental techniques have been invented in order to investigate 
the detailed nature of this non-exponentiality, among them a reduced four-dimensional 
NMR technique [5] , an optical deep bleach experiment [3] and also nonresonant holeburning 
(NHB)|T. Common to all the techniques is that they allow to specifically select a (slow) 
sub-ensemble and afterwards monitor its relaxation. This way the existence of dynamic 
heterogeneities could be verified experimentally In the present paper we adopt the 
definition given in ref.jH] according to which a response or relaxation function is termed 
dynamic heterogeneous if it is possible to specifically address effectively slow, intermediate 
or fast contributions to the ensemble averaged function. Otherwise, the relaxation will be 
termed dynamic homogeneous. It is important to mention that neither this definition nor 
any of the quoted experiments allow to obtain information about any spatial aspects of 
these heterogeneities. 
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NHB consists of a pump-wait-probe field sequence, starting with the apphcation of one 
(or more) cycles of a large amplitude ac field with frequency Q to the sample in equilibrium. 
After a waiting time has elapsed, a small field is turned on to monitor the modified 
response, cf. Fig.l. A qualitative interpretation of the experimental results relies on the fact 
that via the application of the ac field the sample absorbs energy of an amount proportional 
to the imaginary part of the susceptibility evaluated at the pump-frequency f2|7]. One thus 
expects that a frequency selective modification of the spectrum should be feasible only if 
the response is given by a heterogeneous superposition of entities relaxing at different rates. 
This view has proven fruitful in the interpretation of experimental results obtained from 
a variety of samples, including supercooled liquids jH |HI, disordered crystalsj^, amorphous 
ion-conductors ^Oj and spin glasses jUj. A theoretical investigation in terms of a response 
theory has been developed for the case of stochastic dipole-reorientations|12j . which has 
shown that the modified response indeed depends on the absorbed energy, thus supporting 
the above picture P^. 

At this point a note of caution is appropriate. NHB differs fundamentally from the 
so-called spectral holeburning experiment known from nonlinear optics^^, which e.g. has 
successfully been used to investigate the dynamics of two-level systems in glasses [13]. In 
this experiment the optical transition of a dye molecule is altered externally, thus giving 
rise to a hole in the remaining absorption spectrum. In contrast, NHB does not work at 
frequencies allowing to monitor electronic transitions and we are concerned with dynamic 
features which are much slower. 

Because the mentioned theoretical investigation on NHBfT? was concerned with slow 
reorientational dynamics, typically on a time scale of /is ■ ■ ■ s, only purely dissipative dy- 
namics has been considered and inertial effects have been neglected completely. Given 
the fact that the primary relaxation of supercooled liquids is of a dynamic heterogeneous 
nature, the question as to which extent the same holds for the fast dynamics in these 
systems naturally arises. In particular, this point may prove of importance for under- 
standing the physical origin of the so-called boson peak, observable e.g. as low-frequency 
oscillations in nonresonant Raman scattering|16j. On the experimental side, the recent 
developments in far-infrared spectroscopy ^Zj, in particular the application of femtosecond 
Terahertz pulses [T^ have opened the way to investigate the nonlinear dielectric response 
in this frequency range in detail. Thus, there is good reason to hope that the realization 
of an experimental protocol like the NHB field sequence should not be out of reach in due 
course pr?]. 

In this paper we present a theoretical analysis of NHB for the case of damped oscillations 
using the well known model of Brownian oscillators (BO) 14J. This model, which often is 
employed in calculations of the nonlinear optical response in condensed phases, allows us 
to treat both, underdamped and overdamped motions in a coherent way. The outline of 
the paper is as follows. In the next section we describe the model and the theoretical 
aspects of our study. Section HI contains the results of the model calculations along with 
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the discussion, including an estimate of the expected magnitude of the nonhnear effects. 
Finally, we close with some concluding remarks in Sect. IV. 

II. NHB for the Brownian Oscillator Model 

In this section we briefly discuss those features of the BO-model that will be relevant for 
our discussion. Additionally, we calculate the linear and the third-order response functions 
relevant for the NHB experiment. As we are primarily concerned with the response of 
amorphous systems, we have to face the problem of treating spatial correlations of the 
coefficients in a normal-mode expansion of the dipole moments. Therefore a discussion of 
the frequency dependence of the so-called light to vibration coupling constant is included. 

A. General aspects and linear response 

We assume that the normal-modes {q} = {qi, q2, ■ ■ ■} of the system under consideration in 
the THz regime can be approximated as independent oscillators. In a molecular system, 
the qi would correspond to internal degrees of freedom like vibrational or torsional modes 
and in a macroscopic system like an amorphous system the qi represent all vibrational 
excitations including the low-frequency phonon-like collective modes. 

In the classical version of the BO-model with Markovian damping (Ohmic friction) |14j 
and equal masses, rrii = m, the dynamics of the qi is described by the Langevin equation: 

m«. + ™7.4. + =™ = r.(i)-^ (1) 

dqi dqi 

Here, ji denotes the damping constant and Ti{t) is a Gaussian stochastic force with zero 
mean and correlation (rj(t)rfc(t')) = 25i^klil3~^fn5{t — t'), (3~^ = ksT. Furthermore, the 
coupling to the field. Hint, is given byjHj: 

nint = - Jdrfi{r,{q})E{r,t). (2) 

The Fokker-Planck equation corresponding to the Langevin equation, eq.(0), is solved using 
the harmonic the unperturbed model and treating anharmonic terms in ^({q'}) 

as well as Hint in perturbation theory similar to the calculation performed in ref.^2]- 
We mention that the assumption of Ohmic friction can easily be relaxed and a dynamic 
homogeneous scenario can be modelled via a frequency dependent damping 7(0;). This case, 
however, is most relevant for overdamped motions and will be treated in a forthcoming 
publication [201 • The solution of eq.lH)) shows that all quantities of physical interest are 
determined by Ai,i = |(7i + 6i) and Ai_2 = |(7i - 5i) with 6i = {•yf - iuff^'^. 
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In order to calculate the response to an external electric field, we expand the dipole 
moment 

/^(r, {q}) = /io(r) + J2 + ^12 + ■■■ 

and neglect higher order terms as well as cross-terms for simplicity. We will focus on the 
dielectric response but the inclusion of a similar expansion of the polarizability, relevant 
for Raman scattering, poses no problem. Note that we have kept the position-dependence 
of the dipole moment, since in amorphous systems correlations are expected to exist which 
extent over some finite volume, depending on the coherence length of the modes g. |21j. 

In the present model, there are different contributions to a non-vanishing third-order 
response |221 • One nonlinear term stems from the quadratic term in eq. (jHI) • Another possible 
source of nonlinearity has its origin in anharmonic contributions to the potential. In our 
calculations, we use a potential of the form V{{q}) = J2i ^(o'l) with 



V{qi) = muj. 



+ 0{qt) (4) 



where we have scaled the anharmonicity strengths by the harmonic frequency uf. Higher- 
order terms in the potential and also a coupling among the qi can easily be incorporated 
albeit this leads to an increasing number of parameters. 

The polarization following the NHB field sequence for a single mode is calculated in 
0{Ep) and 0{Em)- Without going into the details of the calculations here, we mention 
that the NHB -response Pj^ ^ ^ (t, tnj, tp) is given as a superposition of the ordinary linear 
response and a modification: 

i^*(^)(t,t^,g = pf)(t) + Apf)(t,t^,g (5) 

Here, P^'^\t) denotes the polarization in 0{Em)- In addition, the nonlinear modification, 
AP[^\t,tw,tp), which is 0{Ep ■ Em), contains various terms, to be discussed later. Note 
that the form of the signal is formally identical to the one found in ref.|12^. Thus the same 
phase-cycling as described in ref. |3] can be applied in order to extract the modification 

(X) 

APj . The superscripts 'X' stand for the different response functions measured. In the 
present paper we will discuss three different functions. A field pulse yields the pulse- 
response (X=P), while application of a constant field gives the step-response (X=S), also 
called the integrated response. Furthermore, we consider the response to an oscillatory 
field, denoted by X=AC. Summarizing, we have: 

X=S : EM{t) = EMe{t) ; X=P : EmI^) = ^a/^W ; X=AC : = ^Me"*"* 

where the fields in eq.© are given by E{r,t) = e^^^E{t). 



4 



Denoting the wave vector of the measuring field by k^f , one has for the polarization in 
the linear response regime: 

P^''\t) = ^Em J dvo J t/ri/i(ro):/x(ri):e^^-^^i?f (6) 

Here, p = N/V is the number density of oscillators and -R- (t) denotes the linear response 
of mode 'i'. The latter are different for different time-dependencies of EM{t)- 
In case of the pulse-response we have: 

i?f)(t) = 1 (e-^«.2*_e-^-i*) (7) 

The integrated response function is trivially related to the pulse-response via R[^\t) = 
JodTRPir) and the response to an ac field of frequency lu reads as: 

In the steady state, characterized by e~'^''i/2* = 0, this expression coincides with the fourier 
transform of Rf\t) according to eq.© up to the factor e"**^*. This is exactly the result 
expected from linear response theory [7]. 

B. Nonlinear modifications 

There are various contributions to the cubic response and therefore also to the nonlinear 
modification of the response in a NHB field sequence. These stem from the higher order 
terms in the expansion of the dipole moment, eq.Q, and the anharmonic terms in the po- 
tential, eq.(@]). In the lowest nonvanishing order, there are three relevant terms. One term 
has its origin in the quartic anharmonicity, G4, along with the first order term in the expan- 
sion of the dipole moment. In the following, this term will be denoted as AP/ (t,tw,tp). 
In addition, there is one term, denoted by ^Pi(^harm)i^y^w,tp), which stems from the sec- 
ond order term in eq.Q and the harmonic contribution to V{{q}). Finally, there exists 
one term which mixes the second order dipole moment term with the one of the cubic 
anharmonicity, which we will call AP^f^J^-^^{t,ty;,tp). The principal features of these terms 
are very similar with quantitative differences only. 
The different contributions are given by: 

AP^''\t,t^,tp) = -^^B.EmEI jdv, y'f/ri/i(ro):/i(ri):e^^«^i x 

y< jdv^j rfr3/i(r2):/i(r3):e^'^-(^^+'-^) Ai?f ) (t, t^, t,) (9) 

Here, ARf^\t,tyj,tp) denotes the intrinsic modification of mode 'i' and kp and k^,/ are the 
wave vectors of the pump-field and the measuring field, respectively. 
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From eq.©, the expression for APiflarm)i^^ ^p) obtained by dividing with (— 164), 
replacing /i(rfc)^, /c = 1, 2 by /i(rfc) and using ARf^larm.)it^ ^p) instead of ARf^\t, t^j, tp). 
In order to get the expression for AP.*'"^^^^(t, t^, tp) one has to multiply AP/'^^(t, t^, tp) in 
eq.© by an overall factor (O3/64). Furthermore, AR] {t,tw,tp) has to be replaced by 

j,(X) , 

i,{mix) * 



^^uliix)it't^'tp) and one of the three /i(rfc)-, k^O, by /i(rfc)-'. 



The expressions for the modification Ai?P'*(t, t^, tp) (and also ARf^l^^^^, AR^j^-^^) 
are somewhat more complex than those for the linear response. Here, we concentrate on 
ARl and only mention that the other expressions are of a similar structure without giving 
them explicitly. 

It turns out that for all three measuring procedures, X=S, P, and AC, the modification 
can be written in the form: 

APf)(t,t^,tp) = (^^^{e-2^-^*-x«,2(^^)'/i''(A.,i,A,,2;t) -e-2""*-X^A(^)'^''( 

-e-('''^+'''^)*-X,i(^)x,2(f^)/(A,i, A,2; t)} (10) 
In this expression, the functions Xi^a{^) are defined by: 



with denoting the number of cycles of the pump-field. The functions h {a, b; t) and 
g'^{a, b; t) are linear combinations of exponential functions of time, given explicitly in the 
appendix. We mention that ARf^^\t,tuj,tp) are real quantities despite the fact that the 
Ai,a, « = 1, 2, are complex in the underdamped case. 

From the expression for the modifications AR\ (t, t^, tp), eq. ()10|) . it is evident that the 
absorbed energy plays an important role, albeit the situation is more complex than in case 
of a completely overdamped motion. It can be shown from eq. pup . using the expressions 
given in the appendix, that the modification vanishes in the limit of small as well as large 
times for 7 > and thus is of a transient nature. During the waiting time the created 
non-equilibrium population relaxes due to the terms e"'^'-"*™. This is because there is no 
extra relaxation mechanism in this mo del P^. 

Eq. (jll|) shows that Xi,a{^) fi— >0 and Xi,a{^) ioT Q—>-oo and therefore 

the modification vanishes in these cases also. These considerations show that the nonlinear 
modification of the model considered in ref. and of the BO model are very similar, 
particularly in the overdamped regime, which will be treated in more detail in the next 
Section. 

Because we mainly are interested in oscillatory motions, we consider the signals in the 
frequency-domain as well. Here, the pulse-response is of particular importance as this 
is the quantity that is directly related to the complex dielectric constant in the linear 
regime. Therefore, in addition to the time-domain signals we will consider the fourier 
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transforms of R\ (t) and Ai?- (t, t^, tp) with respect to the measuring time t. However, it 
is important to point out that the hnear relation between the pulse-response and the step- 
response, Rf\t) = JqcItRP (t), does not hold hi the modifications Ai?P^(t, t^;, tp) and 
ARP{t, t^, tp) due to the nonlinear nature of these functions. Similarly, there is no simple 
relation between the fourier transform of AR^ {t, t^, tp), denoted as AR^ {uj, t^, tp) in the 
following and ARl^'^\t,tu],tp). This can be inferred from the definitions of the functions 
h'^'^ and g'^'" , eq. ()A.3|) . For long times the modification measured in the frequency domain, 
AR[^'"\t,tu],tp), vanishes and there is no stationary state evolving as it is the case in the 
linear response, cf. eq.(jHl). This can be seen from eq. ()A.3|) as there is no term of the form 
e~*'^*, but u occurs solely in combinations with Aj^Q. Of course, this behavior also is to be 

expected intuitively because the modification is to be viewed as a transient effect. It is for 

(p) 

this reason that we consider ARj in the frequency domain. 

C. Overdamped limit: Ornstein-Uhlenbeck process 

As discussed above, the relation dP^^\t)/dt = P^^\t) does not hold for the non-linear 
responses. This fact can be shown to hold true for non-linear response functions in gen- 
eral. Here, we will demonstrate it explicitly for the limiting case of a strongly overdamped 
motion, the so-called Ornstein-Uhlenbeck (OU)-process|2n]- All expressions for the OU- 
process can be obtained from the corresponding ones for the BO-model from the lowest- 
order term of an expansion in l/7j. This expansion can either be performed on the Fokker- 
Planck equation itself or on any of the resulting expressions for the linear or nonlinear 
response. For the OU-process, the effective relaxation rate is given by Aj = cu^^ /7j [HI . Ad- 

(S) (S) 

ditionally, the expressions for Rl and Ai?) obtained that way allow a direct comparison 
with those obtained for the model of stochastic dipole reorientations in ref. The linear 
response given in eq.© simplifies to 

RP{t)^ov) = ^e-^^' and Rf\t)^ov) = ^ (l - e~^^') 

For the modifications, however, one finds from eq. lfTIH) : 

ART\t,t^,tp)^ov) = ^,Ain)e-'^^'- (l - e"^^* 

ARf\t,t^,tp)^ov) = ^Ai(n)e-2^»*» (l - e-^^*)'e-^'* (12) 

Here, the pump-frequency dependent amplitude is given by Ai{Q) = [ef (i7)]^ (l — e"^'*''^ 

with ei{uj) = Aj/(Aj — iu). It is evident from these expressions that ARf'\t,tuj,tp) ^ 
jQdrARf \T,ty^,tp). Furthermore, the expression for ARf\t,t^,tp) 

(ou) is similar to the 

one found in ref . |T2| for stochastic dipole reorientations apart from some minor differences, 
which will be discussed in detail elsewhere |20]. 
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If we consider the fourier transform of K- {t)(o\j), {lu)(^o\j) = ^ej(i^) and compare 
it to the signal measured in the frequency domain via apphcation of a field Emg'"^^^, 

for t:^A^^ one finds the standard result of linear response theory, namely -Rj^*"''(t)(ou) = 
e~*'^*i?j-^''(u;)(ou)- However, this does not hold for the modification A_Rj^'"'*(t, t^, tp)(ou)- 
According to eq.()10p. we have: 

which has to be compared to the imaginary part of 

ARf\u,t^,tp)^ov) = -^A(^^)e-'^-*- [Seiico) - e,(^/3)] (13) 

Whereas the dependence on the burn frequency Q is identical, this does obviously not hold 
with respect to the frequency u. The modification AR\^^\t,tuJ,tp)'|^Q■^J•^ clearly exhibits 
a transient behavior with respect to the time t as it vanishes in the limit of short and 
long times. There is no stationary state which would allow to set a time-window for 
observation, in marked contrast to the linear response. These considerations show that it 
is most meaningful to consider ArP{uj, t^, tp), i.e. the fourier transform of ARP{t, t^, tp), 
in the frequency domain and not the signal measured via application of an ac-field. 



J sin (iut) 



D. Light to vibration coupling 

As already noted above, in amorphous systems spatial correlations of the dipole moments 
are expected to exist which extent over some finite volume, depending on the coherence 
length of the modes g^. Therefore, the position-dependence of the /u(r). and the /u(r)f 
plays an important role. This fact complicates the relation between the intrinsic response 
functions and the respective polarizations, eqns.® and (0). Thus, the situation is very 
similar to the case of Raman scattering in glasses|21^. 

In the expression for the linear response, eq.®, we have to consider the correlation 
function Jdvo JdTifi{ro)[fi{ri)[, where we have used e^^'^'^^ ~ 1, i.e. the k — > limit of 

k)^/x(k) ■). In order to proceed in the calculation of this function, we make the same 
assumptions as they are typically used in calculations of the Raman-scattering intensity 
from amorphous systems j2Il 1211 • assume that the coherence length of the mode qi is 
much smaller than the wavelength of the light. Additionally, we disregard the dependence 
of all quantities on the polarization of the mode and of the light. In an amorphous system 

k) ■/i(k)^) is expected to have a broad flat maximum in the vicinity of k = 0, in vast 
contrast to the situation in crystals. This is the reason for the breaking of the momentum 
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transfer selection rules in the case of Raman-scattering[21j. A comparison of our expression 
to the ones obtained for the Raman-scattering intensity shows that k)^/i(k)^) is pro- 
portional to the so-called light to vibration coupling C{uJi). Indeed, a very similar behavior 
of the Raman- and Infrared-coupling constants has been observed in a silica glass j2Sl- We 
thus assume k)^/i(k) •) r^C^Ui) apart from an overall prefactor which replaces 

la^P occuring in the Raman case, i.e. 

{f,{-kyxm^\f^Tc{u:.) (14) 

In the expression for the modification AP/^\ eq.Q, the fourier transform of the four- 
point correlation function (/i(ro)^/i(ri)^/i(r2)^/i(r3)^) occurs. In general such quantities are 
extremely difficult to calculate. Therefore, in order to be able to proceed we assume 
that the probability generating functional for the correlations can be represented by a 
Gaussian. Note that this is exact for harmonic vibrations and otherwise represents a mean- 
field approximation for the spatial correlations of the modes g^. With this approximation 
the four-point correlations factor into products of two-point correlations yielding 

Jdro Jdr, jdv2 Jdr,f,{royMrMr2%f^ir3%e'^''''^^^^^ 

in the expression for APj^\ eq.®- 

For API$'_, another four-point correlation function, (Mr„)MrO;V(r,)'.>(r,);) oc- 
curs, cf. the discussion in the context of eq.Q- In the mean-field approximation one 
thus has to deal not only with k)^/i(k)^), but also with the unkown correlation func- 
tions (yu(— k) •/i(k) •') and k) •'yu(k)"). For the function k)^/i(k)") it is reasonable 
to assume that it is vanishingly small in isotropic systems because of the odd number of 
derivatives. We thus assume 

(^(-k):^(k):') ^ (15) 

For the function k)^'/i(k)^') one expects a similar behavior as for k) •/i(k) •). In 
general, however, the frequency-dependences may be different. Therefore, we write 

(/.(-k):v(k):')^i^'/pc'(c.,) (16) 

and this way find 

(/i(ro):/^(ri):V(r2):V(r3):) ^ C(.;,)C"(.;,) 

Only if one assumes in an ad hoc manner that the /i(r)-' couple to the same elasto-optical 
constants as the /i(r)- one has C'{uJi) = C{uJi) in the expression for AP- j-^^^^^. Due to the 

(X) 

lack of a theoretical argument in favor of such an approximation we mainly focus on AP{ 
throughout the present paper. 

The last term occuring in eq.© is the 'cross-term' A.Pl^^-^y The relevant four-point 
correlation function in this context is (yu(ro)i/u(ri)^yu(r2)^';u(r3)i)- However, in every term 



9 



of the mean-field approximation the two-point correlation function k) ■yu(k)") occurs. 
Therefore, due to eq. (fTH|) we approximately have 



Various models for the low-frequency excitations in glasses have been used to calculate 
C{uJi). All of them yield power laws, C{uJi) ~ cj", with exponents n ranging from ?t, = in 
the soft potential model[2n] to n = 2 in the harmonic modelj21]. Furthermore, when C{uJi) is 
analyzed via a comparison between experimental Raman- and neutron-scattering intensities 
from the same sample |i27j often a C{uJi) ~ Ui dependence is found. Also a nonvanishing 
limiting value C(a;j — *• 0) 7^ has been reported|2Sj. It also should be mentioned that 
even an explicit k-dependence of the low-frequency Raman-scattering intensity has been 
observed in a silica glass|2n]- Thus, given the uncertainty regarding the light to vibration 
coupling, in the present paper we will use C{u!i) ~ cj" with n = 0, 1, 2 in order to discuss 
the possible behavior. 

In actual calculations, we furthermore disregard the dependence of /i(r)- on the mode- 
index i and write for the coupling constant 



It should be noted that if the system under consideration consists of independent particles, 
spatial correlations of /i vanish and the relation between the polarization and the response 
functions become trivial. In this case we simply have C{u!i) = C'{uJi) = 1 and only the 
corresponding prefactors and occur in the expressions for the polarization. 

E. Polarization 

In addition to the approximations discussed above we assume that the parameters of the 
BO-model, 7, and Ui, are distributed according to some distribution functions g{uji) and 
piji), which we choose to be independent. For the density of states (DOS) of the low- 
frequency vibrations we will use g{u!i) ~ u>^, m = 2,4 as they follow for the Debye model 
and the soft potential model, respectively. For computational convenience we additionally 
introduce a high-frequency cut-off Uc. 




(t, tyj, tp) — 



{fi{-k)[fi{k)[) = \fifC{u;,) with Ciu,)ocu', 



(17) 




(18) 



With these approximations and definitions we have from eq. © : 



m J J 



(19) 



and an analogous expression for the modification, eq.©: 




(20) 
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As already noted in connection with eq.Q and in the previous section, the expres- 
sion for APlj2rm) obtained from the one for AP^^^ by replacing [— |04|/i'|'^C(co'j)^] by 
[\lJ,'\'^\lJ,"\'^C{(jJi)C'{(jJi)] and using ARf^l^^^^ instead of AR^^K These expressions will be 
used to discuss the response of a collection of Brownian oscillators to the NHB pulse se- 
quence. 

At this point it should be pointed out that in case of G4 > the sign of the two functions 
AP(^) and AP'^^^^ ^ will be the opposite. This is an example of the fact that one cannot 
predict the sign of a third-order response in general. 

III. Results and Discussion 

In this section we will discuss the results of model calculations and show that indeed 
the NHB field sequence is capable to detect dynamic heterogeneities provided they exist. 
In the discussion of the nonlinear modification we will mainly focus on AP^^^ but also 
show results for APjj^^^^y In addition we will give an order of magnitude estimate of the 
expected effects, as this is of utmost importance for possible experimental realizations. 

A. Dielectric loss 

We start with a brief discussion of the dielectric loss, t"{uj). From eq. ()19j) and the fourier 
transform of eq.© one obtains e"{uj) = {eoEM)~^y^[P^^\t)], where eo is the permittivity 
of free space. For simplicity we restrict ourselves to a single damping constant, i.e. we use 
p{li) = — 7) in eq. (fT9|) . We then explicitly have: 

= ^l/i? / dujM^,)C{u;,)—^ ,2 (21) 

The limiting behavior of €"{uj) can easily be extracted from this expression. For vanishing 
7 eq. (PT|) reduces to e"{uj) = p{eom)~^ {it /2)\fi'\'^C (u) g{u) /u which is directly proportional 
to the Raman intensity scaled to the Bose factor. Using eq. (fT7|) for C{uJi) and eq. (fTH|) for 
g{uJi) this implies e"{uj) oc u™'^"'^^. For finite 7 we have 7 > at low frequencies and 
therefore ^"{uj) ~ uj. For frequencies u; > 7 this behavior changes into that of undamped 
oscillations. Since the DOS has a high frequency cut-off ujc we find for u > Uc e"(^) ~ 
in the overdamped case and €"{uj) ~ uj~^ otherwise. This behavior is exemplified in Fig. 2, 
where we plotted e"{uj) versus u /uc for m = 2, n = 1 and various values of the damping 
constant 7. Also included in Fig. 2 as the thin dotted line is the dielectric loss for a Debye 
relaxation, e"j^{uj) = ■yu/^u^ + 7^), which behaves as e^(c<j) ~ •j'^u and e'^iuj) ~ '^uj~'^ 
for small and large cj, respectively. This shows that in the overdamped case we have an 
apparent distribution of relaxation times, giving rise to a sublinear increase of t"{uj) for 
small UJ. Of course, the cusp in e"{uj) for small 7 has its origin solely in the assumed 
high-frequency cut-off of g{uj). 
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B. Nonlinear modifications 

Before we discuss the behavior of the modification AP^^\uj,tw,tp) it is instructive to 
consider the response associated with a single mode g,. In Fig. 3 we show the real and 
imaginary parts of AR^^^ {^,tu,, tp), the fourier transform of Ai?-^^ (t, t^, tp) given in eq-lfTIH). 
for a burn frequency fl = 1.0. The full lines correspond to an underdamped oscillator 
whereas the dashed line represents the overdamped case. Note that Xi^a{^) in this situation 
only determines the overall amplitude of the modification. Only when a distribution of 
modes is considered the frequency selectivity of Xi,a{^) becomes important. In the present 
situation of a single mode this only means that the overall amplitude of ARf\t,tw,tp) 
is changed if the burn frequency Q is varied. In the underdamped case a clear resonance 
occuring at uj = Ui is visible in both, A_R-^^(u;, t^, tp)" and ARf\uj,tw,tpy. Additionally, 
a weak resonant behavior is found at a; = Scjj. This can be understood from the definition 
of and g^, cf. eq. flA.ljl . and stems from terms of the form l/[3Aj_Q, — ioj]. On the other 
hand, in the overdamped case the maximum modification is found at u; ~ Aj (more precisely 
uj ~ 0.86Aj). In the limit of the OU-process, the result for ARf^ has been given above in 
eq. (fT3|l . The factor uf in the denominator explains the smallness of ARf\uj,tuj,tp)[ov), 
which in Fig. 3 was multiplied by a factor 10*. In this case of overdamped dynamics the 
imaginary part appears to be better suited to find the maximum modification. Thus, in 
the following we concentrate on ARP{u,tu„tp)" although in some situations the real part 
may also yield helpful information. 

The behavior of ARf\uj,ty^,tp) as a function of uj also determines the modification 
of the polarization, AP^^\iu,t^,tp) according to eq. lj^Uj) . In the upper panel of Fig. 4a we 
plotted AP^^\uj)" = AP^^\uj, t^,, tp)" versus frequency for various pump- frequencies and 
one cycle of the pump-field. As in Fig. 2, we chose m = 2 and n = l, i.e. g{uji) ocuf, C{LJi) = 
uji. Additionally, we consider a single damping constant, thus writing p(7j) = 5(7j —7) with 
7 = 10-2tuc in eq.(I2ni)- We find that for all VL the minimum of the modification is observed 
at frequencies tUmin ~ In the lower panel AP'^l^^-^{uj)" is shown. In this calculation 
we assumed C'ioJi) =C{uJi) for simplicity and used the same parameters as in the upper 
panel. It is evident that the main differences are the reversed sign of the modification and 
that Wmax is slightly larger for all pump frequencies. Additionally, in the underdamped 
regime a slight minimum is observable at C(j = 7, which is absent in AP^^\uj)" . However, 
the most important conclusion that can be drawn from Fig. 4a is that the modification of 
the response clearly exhibits a pronounced f2- selectivity. This means, that NHB is capable 
to 'detect' dynamic heterogeneities also in the regime of oscillatory motions. 

In addition to the f2-selectivity we can discriminate between underdamped and over- 
damped modes due to the functional form of AP^^\ijj)" . In the overdamped regime the 
modification only shows a minimum, cf. eq. (fT!?|l (a maximum in case of AP'^l^^ )('^)")- 
the underdamped case an additional dispersive behavior is observed. This means that also 
in the macroscopic polarization the main features already observed for the response of a 
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single mode persist. In Fig.4a, the dashed hnes correspond to overdamped modes because 
here <^ 7 and thus mainly the response of those modes with Aj ~ f2 is modified by the 
pump. For higher burn frequencies, however, mainly the response of underdamped modes 
is modified, cf. the full lines in Fig. 4a. Here, the maximum modification is around uj^^uJi 
and therefore ci;j~fi. Of course, for other choices of the damping constant 7 the cross-over 
from underdamped to overdamped modes takes place at a different frequency. 

If we compare the form of the modifications in the underdamped regime, the full lines 
in the upper panel of Fig.4a, to those for a single underdamped mode shown in Fig. 3, it is 
evident immediately that the range of non-vanishing modification extents over a consider- 
ably larger tu-range. This has its origin mainly in the functional form of cf. eq. (jll|) . 
in particular the factor (1 — e~'*'» °*p), which shows oscillatory behavior in the underdamped 
case. At low frequencies some more strongly damped modes contribute to the signal, ren- 
dering the modification broader as compared to Fig.3. This effect can be suppressed by 
increasing the number of cycles of the sinusoidal pump-field. The condition e"^' "**" <^ 1 
is fulfilled for A^^i7/(7r7j) and thus for a large number of cycles the behavior of a single 
oscillator is approached. Therefore, in the regime of underdamped modes the number of 
cycles A^ may turn out to be a useful experimental parameter. This is in contrast to the 
situation of overdamped modes where increasing A^ mainly gives rise to a small increase in 
the overall amplitude due to the factor (1 — e~^'*p). 

Next, we consider the frequency at the minimum modification, d^min. This quantity 
is plotted versus VL in Fig.4b for two pairs of {m,n). Remember that these parameters 
determine the effective distribution of uJi contributing to the signal because according to 
eq.(|^ we have (yf(ci;j)C(ciJj)^ occu™"'"^". Due to this behavior, the equivalences (2, 1) = (4, 0) 
and (2, 2) = (4, 1) hold for the pairs (m, n) shown. As in Fig.4a, the upper panel shows 
ciJinin extracted from AP(^)(ti;)" and in the lower panel ujmax obtained from API^^^^^-^{uj)" 
is plotted. In the latter case, we again assumed C'{ui) = C{LJi). It is evident from Fig.4b 
that cUminOcfi. Only in the regime of overdamped motion different pairs (m,n) give rise to 
some minor differences in the values of u^in (t^^max)- Of course, the hnear relation cUminOcf] 
has its origin in the form of the effective distribution of Ui in the oscillatory regime and the 
corresponding distribution of Aj = ujf/'~f in the overdamped regime. In the underdamped 
case, one does not only have a minimum in AP^^\uj)" but also a maximum in addition, 
cf. Fig. 4a. However, as discussed above, with increasing number of cycles, the differences 
between the minimum- and maximum-positions diminishes. 

Sofar, we have considered the case of a fixed damping constant, which of course should 
be interpreted as an average value. It is evident from the discussion of Fig.4 that the 
value of 7 can hardly be extracted from AP^^\uj,tw,tp)" and this does not change if a 
distribution of damping constants is assumed in the calculations. On the other hand, it 
is known from previous studies that AP^^\t,tu,,tp) is sensitive to relaxation rates and we 
thus expect the same for AP^^\t,tw,tp). In order to demonstrate that in the oscillatory 
regime measuring AP^^\t^t^^tp) may provide information about the damping constants. 
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in the upper panel of Fig. 5 we plot this quantity versus measuring time for the same 
parameters as used in Fig. 4. Here, we used VL/uJc = 0.1 and various numbers of cycles of 
the pump-field, N. The first thing to notice is that for = 1 hardly any oscillations are 
visible in the time-domain, whereas this behavior changes with increasing A^. This fact has 
the same origin as the decreasing broadening of the modifications for larger N, as discussed 
in context with Fig. 4a. For a single cycle of the pump-field AP(^)(co', t^, tp)" shows a very 
broad resonance and accordingly the dot-dashed line shows almost no oscillations. For 
= 10 AP^^\t,tw,tp) shows the expected oscillatory behavior. More important for the 
present discussion, however, is the fact that for a large number of cycles it is evident that 
the signal has its maximum envelope around t~7~^. This is demonstrated by the dotted 
lines, which are of the form txe""'^. This functional form of the envelope can be derived 
analytically from the expression for the pulse- response of a single mode, eq. ffTIUl . 

The lower panel of Fig. 5 shows AP^^^j^,^ -^(t, t^, tp) for the same parameters and the 
additional assumption C'{uji) = C{uJi). In this case, however, the situation is quite different. 
Here the main effects are observed at 'yt = Q and the oscillatory behavior around 7t ~ 1 is 
strongly suppressed. On the other hand, the maxima of the signal are by about a factor 
of ten smaller than in case of AP^^\t,tw,tp). (The envelope function (dotted line) is the 
same as in the upper panel but divided by a factor of 30.) Thus, we conclude that also in 
the oscillatory regime NHB experiments may be useful in order to get information about 
the damping, although this might require to perform experiments with various numbers 
of pump-field cycles. In particular, this should be feasible if the main contribution to the 
signal stems from AP^^\t,t^,tp), i.e. for very small |yu"p. 

C. Estimated magnitude of nonlinear effects 

For an experimental realization of the NHB experiment in the THz regime it is of course 
desirable to have at least an order of magnitude estimate of the expected nonlinear effects. 
For this reason we consider the ratio AP/P, where P and AP are shorthand notations for 
the linear response and the nonlinear modification. Experimentally, the sum P* = P + AP 
will be observed and AP has to be extracted by a proper phase cycling procedure, cf. 
ref.Hj. 

We start our consideration with AP{u!)" for a single overdamped mode, given in eq.()13p. 
For our rough estimate we consider the frequency of the maximum modification, i.e. fl = 
Ai = u, which allows us to neglect e-'(Aj/3) as compared to e-'(Aj) and to use exp(— Ajtp) ~0. 
Furthermore we set the waiting time to zero. This way we find: 

AP, _ 3E^U4l 

D — Q 2 4 4 

Pi 8 m^tUj 

In order to proceed we replace the derivative of the dipole moment by fi[ ~ {6^/6qi) and 
assume that {6qi)'^ is either related to the thermal expectation value (qf) = [ksT) / {muf) 
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or the average displacement q^. This allows us to write 



AP 3 (5^-E_ 



2 

P \ ^2 



P 8 \ kuT 



qfe, (22) 



In order to proceed, we need an estimate of the quartic anharmonicity B4. For this purpose 
we utilize the detailed comparison of the predictions of the soft-potential model (where a 
potential of a similar form as in eq.(j3]) is used) with experimental results on a variety 
of different glassforming liquids [30] ■ From this comparison we find an average value of 
04 ~ 2 ■ lO^^m"^. Assuming a value of on the order of IDb/lAiSl], the value of the 
actual displacement and the one to be chosen for 6fi are not independent of each other. If 
qi = xx lA we should use (5/i~xx IDb. Note that this gives an estimate for g,^G4~ 2x^-10^, 
which should be smaller than unity in order for the perturbation expansion used in Sect. II 
to be meaningful. Assuming a pump-field amplitude Ep = z x lO'^kV/cm one then finds, 
using (10^kV/cmxlDb)/(A;ijxlOOK)~ 2.5, that AP/P~ Thus, for a reasonable 

value of x~0.05 one has AP/P~ 3 ■ 10^^ z'^, meaning that AP/P is in the percent range 
for z~2 or Ep~2 • lO^kV/cm. 

In this context, it should be mentioned, that the above estimate of course should not 
be used if one is concerned with the primary relaxation of high temperature liquids. This 
is because in that case the dielectric relaxation is mainly determined by the mean square 
fluctuations, (A/x^), of the static molecular dipole moment /i due to the molecular tumbling 
motion. In our formulation this means that we loosely should identify (Sfiqi) G4 with 
(A/x^), cf. ref. ^21- Assuming isotropic rotational jumps as a simple but realistic model for 
the stochastic reorient at ional motion in liquids [^2] allows us to identify qi with a 'rotational 
jump length', q^ ~ Rh sin {(prot) ■ Here, the hydrodynamic radius Rh is on the order of 
lA[HH!, yielding qfQi ~ 6, if the mean jump angle (prot is on the order of ten degrees. 
Again we write Ep = z x lO^kV/cm in eq. ()22|) . but now we use T = 300K. For a dipole 
moment of IDb this gives (AP/P)^^.^^ ggg^,^ ~4z^. Therefore, a pump-field amplitude of 
Ep ~ 500kV/cm is sufficient to get an effect in the percent range. This value is in harmony 
with the amplitudes used in the NHB experiments on supercooled liquidsjl]. 

Next, we consider the more relevant case of underdamped oscillatory motion. In order 
to give an estimate of the ratio APiuj)" / P{uj)" in this regime we now consider a distribution 
of eigenfrequencies uoi similar to the discussion in context with Fig. 4a (upper panel) and 
C{ijJi)o:.uJi. The expression for the dielectric loss, eq. ipTj) . allows us to write 



P{uj)" 



tot {uj)p^^ EpQ4- 7- 72 (23) 

Jduig{uJi)uJiRl \uj)" 



In this expression, pm = 'mpc:^10^kg/m^ is the samples mass density. In order to obtain a 
peak in the THz-range we now use a DOS of the form g{uji) (xujf exp(— rtUj), with r = 4/3 ■ 
10""^^ Hz~^, yielding a peak at 0^1 = 1.5THz. Consequently, we set VL= 1.5THz. Additionally, 
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we fix the damping constant to value of 7 = 0.1THz. We point out, tliat a variation of 7 
lias only a minor effect on the resulting amplitudes as long as the motion is in the strongly 
underdamped regime. For the dielectric loss in the THz-range we approximately have 
e"(THz) ~1, as found for glyceroljM]. Evaluating the integrals in eq.(j2SI) at c^min — ITHz 
(cf. Fig. 4) allows us to calculate AP^ax- For 64 we use the same value as before, 64 
lO^^m^^ and find 

AD TT-S 

3.10-i5_^ 



max 



P 



(ITHz) 



This way we obtain APmax/-P — 0.3 for a pump-field amplitude of Ep = lO^kV/cm. These 
considerations allow us to conclude that - with the assumptions used - the ratio AP/P 
should be larger in the oscillatory regime than in the overdamped regime. Of course, from 
an intuitive point of view this does not come as a surprise because some resonant behavior 
even is expected in this case. 

Sofar, we have considered the nonlinear modification AP{uj)". We close this section 
with a brief discussion of the relative magnitudes of AP and AP(^harm.), the nonlinear 
modification stemming from the second order term in the expansion of the dipole moment, 
eq.(jni). As is evident from Fig. 4a, the magnitudes of the normalized effects are quite 
similar. Thus, we have to compare the prefactors 304|/i'p and For this purpose we 

proceed in exactly the same way as we did in the discussion of the overdamped motion. 
We approximate by and by ((5[(5/x])^/(3(g^)g^). This way we find for 

the ratio 

If we now assume that 6fi^x x IDb and simply by analogy that S[6fi] x IDb, we have 
364|/i'p~2 ■ 10^|/i"p and thus a negligible effect due to AP(^harm.)- However, it should be 
borne in mind that apart from the necessity of estimating 5[6fi] also the function C'{uji) 
occuring in the expression for AP(^harm.) is unknown. On the other hand, an experimental 
determination of the relative relevance of AP and AP(^harm.) appears feasible due to the 
fact that the effects are of opposite sign. 



IV. Conclusions 

We have calculated the response of a system consisting of a collection of Brownian oscilla- 
tors to the NHB field sequence. For simplicity, we have restricted ourselves to a classical 
calculation and assumed Ohmic friction, i.e. time-independent damping constants. A more 
general case will be treated in a separate publication [201 • Furthermore, we have restricted 
the calculations to a quadratic expansion of the dipole moment in the normal-mode co- 
ordinates Qi and included no higher terms than the quartic anharmonicity (oc G4) in the 
potential. We did not allow for any coupling among the modes. We do not expect that 
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any of these assumptions presents a severe restriction on the apphcabihty of our results. 
Additionally, we expect a similar behavior in case that the light couples predominantly to 
the polarizability instead of the dipole moment. 

We have applied our calculation to the situation encountered in amorphous systems 
in the frequency range of the so-called boson peak. We do not claim that a BO-model 
quantitatively describes the features of the vibrational excitations in this frequency regime. 
However, there appears to be some consensus about the quasi-harmonic nature of these 
vibrations. Although different models for the dynamics in the boson peak regime attribute 
different damping mechanisms |2 6 1 EHl ISE]; experimental results often can be fitted to a 
simple damped oscillator susceptibility function j^Tj. Because less is known about the 
detailed form of the DOS in the THz range, we simply parameterized it as g{u}) oc uj'"^. 
In contrast to the situation of low-frequency Raman scattering, we are concerned with 
a nonlinear response and therefore have to face four-point spatial correlation functions 
of some elasto-optical constants. In a simple mean-field treatment we factorized these 
four-point correlation functions. For the various contributions to the nonlinear response 
different two-point correlation functions are relevant. We have focussed predominantly on 
the contribution that depends on the square of the so-called light to vibration coupling. For 
the latter, we used a functional form C(a;) =uj^ with various values for n in order to be as 
flexible as possible and not to rely on some specific model. Along with the behavior of the 
DOS we thus can investigate the behavior for models ranging from Debye-like phonons|24j 
to soft-potential quasi-harmonic modesl26]. 

Our main finding is that also in the regime of oscillatory motions is NHB able to 
'detect' dynamic heterogeneities via a pump-frequency selective modification of the linear 
response. In contrast to slow relaxation processes, where usually the response to a step 
field is monitored in the time domain it seems to be advantageous to consider the fourier- 
transform of the response to a field pulse in case of underdamped motion. For all relevant 
pairs (m,n) we find a linear relation between the frequency of the maximum modification 
and the pump-frequency Q. For underdamped modes oscillators with uJir^Q are addressed 
primarily whereas in the overdamped case those with Aj~fi are most relevant. The latter 
situation is similar to what has been found in earlier investigations. If the dynamic range of 
the experiment is large enough, a cross-over from underdamped oscillations to overdamped 
relaxation should be observable in principle. 

An important point in the context of nonlinear response functions regards the relation 
between the pulse- and the step-response. Whereas these functions are trivially related in 
the linear regime this has been shown not to be the case for the nonlinear modifications. 
We explicitly demonstrated this behavior for the strongly overdamped case, i.e. the OU 
limit. For this case a modification of the step-response has been found which is very similar 
to the expression given earlier for a model of stochastic dipole reorientations 

Regarding the magnitude of the expected nonlinear effects, we estimate that pump-field 
amplitudes on the order of 10^ ■ ■ ■ lO^kV/cm should suffice to find nonlinear modifications 
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which are on the order of one percent of the hnear response of the system. With the newest 
technical achievements, this leads us to conclude that the experimental realization of NHB 
and other nonlinear experiments in the THz regime should be practible. 

In conclusion, we have shown that NHB experiments in the THz region are expected 
to yield valuable information about the physical origin of the long-discussed vibrations 
around the boson peak. In particular the question as to which extent the dynamics in this 
regime has to be viewed as dynamic heterogeneous can be answered, thus shedding further 
light on the nature of the dynamics in disordered systems. 
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Appendix 



In this Appendix we give the exphcit expressions for the modification of the response, 
ARf^\t,ty^,tp), that are used in the text. 
Using the abbreviation 



a — b 

the functions h^{a, b; t) and g^{a, b; t) are given by: 

h''{a,b;t) = l[^3,.„(t)-^,+26;6(t)] 



gP^a, b; t) = [U2b;a{t) ' 6+2a;6(t)] (A.l) 

a + 



and 



h^{a, b;t) = — [6fe;a(i) + ^a+26;2fe(t) " ^3fe;a(t) " ^a+2b;bi*)] 



g\a,b-t) 



a + b 



h^^{a,b;t) = ^ 



^a+b/2;b/2{'^) ^6+a/2;a/2 1'^) 



i^2b (t) — ^2b+iuj,b(t)) 



a + b 



3b — a 



{^2b+iuj,a(t) — C,2b+iuj,3bit)) 



9^''{a,b;t) 



a + b 



+b+iuj,2a+b 



it) - e 



a+b+iLxj,b 



it)) 



--^ iCa+b+iuj,2b+ait) — ^a+b+iuj,ait)) 



(A.2) 



(A.3) 



In order to calculate the fourier transform of -R-^^(t) and ARf\t,tw,tp) it is sufficient 
to note that these can easily be obtained from eqns. ^ and (fTUI) by noting that 



a — b \a — iuo b — iu 
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Figure captions 



Fig.l : The field sequence for the nonresonant hole burning (NHB) experiment: One or 
more cycles of a strong sinusoidal field Ep{t) = Ep sin (Qt) are applied to a sample 
in thermal equilibrium. After a waiting time ty^ the response to a small field EM{t) 
is monitored. 

Fig. 2 : The dielectric loss €"{uj) according to eq. ()21|) versus frequency for g{uji) = 
{?)/uj1)ujf6{(jJc — uji), C{uJi) = oji and p/(meo)|/i'P = 1. The damping constant 7 
is chosen as 7 = 0,0.01,0.1, 1.0, 100 and a;c = 10. from bottom to top. Also shown is 
e'lj{uj) for a Debye relaxation (thin dotted line) for comparison. Note that for 7 = 100 
the motion is strongly overdamped in the whole frequency range. 

Fig. 3 : Imaginary (upper panel) and real (lower panel) part of ARf\uj,tw,tp) (abbre- 
viated as ArP{uj)) versus measuring frequency u for Ui = 1.0, 7^ = 0.2 (full lines; 
underdamped case) and cUj = 10.0, 7^ = 100.0 (dashed lines; overdamped case). One 
cycle of the sinusoidal pump-field with frequency f2 = 1.0 has been used and the wait- 
ing was set to zero, t^^ = 0. ARf\uj, t^, tp) is obtained from Ai?|^^(t, t^, tp), eq. ipTH) . 
via fourier transform as explained in the text. The results for the overdamped case 
are multiplied by a factor 10*^. 

Fig. 4 : a: AP^^\uj)" versus to for various burn frequencies f2 (upper panel). The 
DOS and the light to vibration coupling constant are the same as in Fig. 2. The 
other parameters are {3p/m^Q4^EMEp\fj,'\^) = 1, = and 7 = lO'^cUc. The 
burn frequencies are chosen as Q/uc = 0.0003,0.001,0.003 (dashed lines) and 
fl/uc = 0.01,0.03,0.1,0.3 (full lines). The dashed lines present overdamped oscil- 
lators while the full lines correspond to underdamped modes. In the lower panel 
'^■P(harm)i^y' plottcd, wherc additionally C'{uji) =C{uJi) is assumed and we set 
{p/m'EMEl\pf\p'f) = 1- 

b: Upper panel: The frequency at which AP'^^\uj)" takes its minimum value (cf. 
Fig. 4a), cjmin/i^c versus burn frequency Q. Lower panel: u^ax/uc versus Q for 
^^{harm)i^)" ■ '^^^ parameters for the DOS and the light to vibration coupling 
constants are written as g{uj) oc cj™" and C{uj) = u^. The results shown are for two 
pairs {m,n) and a damping constant 7 = lO^^cJc- 

Fig. 5 : Ai?(^^(t, t^, tp) (upper panel) and AR^^l^^-^{t,tw,tp) (lower panel) versus mea- 
suring time for the same parameters as used in Fig. 4, Q/uJc = 0.1 and various numbers 
of pump-field cycles, A^. The dotted envelopes are of the form te~'^*. In the lower 
panel, the envelope has been divided by a factor of 30. 
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Haberle and Diezemann; Fig.1 
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Haberle and Diezemann; Fig. 2 
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Haberle and Diezemann, Fig.3 
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Haberle and Diezemann, Fig. 4a 
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Haberle and Diezemann, Fig. 4b 



28 




Haberle and Diezemann, Fig. 5 
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